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We theoretically investigate the deflection-induced coupling of an electron spin to vibrational 
motion due to spin-orbit coupling in suspended carbon nanotube quantum dots. Our estimates 
indicate that, with current capabilities, a quantum dot with an odd number of electrons can serve 
as a realization of the Jaynes-Cummings model of quantum electrodynamics in the strong-coupling 
regime. A quantized flexural mode of the suspended tube plays the role of the optical mode and 
we identify two distinct two-level subspaces, at small and large magnetic field, which can be used 
as qubits in this setup. The strong intrinsic spin-mechanical coupling allows for detection, as well 
as manipulation of the spin qubit, and may yield enhanced performance of nanotubes in sensing 
applications. 



Recent experiments in nanomechanics have reached 
the ultimate quantum limit by cooling a nanomechani- 
cal system close to its ground state jl . Among the vari- 
ety of available nanomechanical systems, nanostructures 
made out of atomically-thin carbon-based materials such 
as graphene and carbon nanotubes (CNTs) stand out due 
to their low masses and high stiffnesses. These properties 
give rise to high oscillation frequencies, potentially en- 
abling near ground-state cooling using conventional cryo- 
genics, and large zero-point motion, which improves the 
ease of detection 2, 3 . 

Recently, a high quality-factor suspended CNT res- 
onator was used to demonstrate strong coupling be- 
tween nanomechanical motion and single-charge tunnel- 
ing through a quantum dot (QD) defined in the CNT 
[3]. Here, we theoretically investigate the coupling of a 
single electron spin to the quantized motion of a dis- 
crete flexural mode of a suspended CNT (see Figjl]), 
and show that the strong-coupling regime of this Jaynes- 
Cummings-type system is within reach. This coupling 
provides means for electrical manipulation of the electron 
spin via microwave irradiation, and leads to strong non- 
linearities in the CNT's mechanical response which may 
potentially be used for enhanced functionality in sensing 
applications [SHZ]. 

In addition to their outstanding mechanical properties, 
carbon-based systems also possess many attractive char- 
acteristics for information processing applications. The 
potential for single electron spins in QDs to serve as the 
elementary qubits for quantum information processing 
[8] is currently being investigated in a variety of systems. 
In many materials, such as GaAs, the hyperfine interac- 
tion between electron and nuclear spins is the primary 
source of electron spin decoherence which limits qubit 
performance (see e.g., [5]). However, carbon-based struc- 
tures can be grown using starting materials isotopically- 
enriched in ^^C, which has no net nuclear spin, thus prac- 
tically eliminating the hyperfine mechanism of decoher- 




FIG. 1: Schematic of a suspended carbon nanotube (CNT) 
containing a quantum dot filled with a single electron spin. 
The spin-orbit coupling in the CNT induces a strong coupling 
between the spin and the quantized mechanical motion of the 
CNT. 

ence [TU], leaving behind only a spin-orbit contribution 
112j . Furthermore, while the phonon continuum in 
bulk materials provides the primary bath enabling spin 
relaxation, the discretized phonon spectrum of a sus- 
pended CNT can be engineered to have an extremely 
low density of states at the qubit (spin) energy splitting. 
Thus very long spin lifetimes are expected off-resonance 
|13j . On the other hand, when the spin splitting is nearly 
resonant with one of the high-Q discrete phonon "cav- 
ity" modes, strong spin-phonon coupling can enable qubit 
control, information transfer, or the preparation of entan- 
gled states. 

The interaction between nanomechanical resonators 
and single spins was recently detected [13] , and has been 
theoretically investigated [HI [TB] for cases where the 
spin-resonator coupling arises from the relative motion 
of the spin and a source of local magnetic field gradients. 
Such coupling is achieved, e.g., using a magnetic tip on a 
vibrating cantilever which can be positioned close to an 
isolated spin fixed to a nonmoving substrate. Creating 
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FIG. 2: Energy levels of the four dimensional (due to spin 
and valley) orbital ground state subspace of the QD, as a 
function of the magnetic field parallel (-By) and perpendic- 
ular (-B_l) to the CNT axis. The boxed areas indicate the 
working regime for the spin qubit (S) and Kramers qubit 
(K), the latter being operated either in a longitudinal (Kz) 
or perpendicular (Kx) magnetic field. Parameter values [30j : 
Aso = 170 ^eV, Akk' = 12.5 /^eV, /Xorb = 330^eV/T. 

strong, well-controlled, local gradients remains challeng- 
ing for such setups. In contrast, as we now describe, in 
CNTs the spin-mechanical coupling is intrinsic, supplied 
by the inherent strong spin-orbit coupling |17H20j which 
was recently discovered by Kuemmeth et al. [21]. 

Consider an electron localized in a suspended CNT 
quantum dot (see Fig. [ij . Below we focus on the case of 
a single electron, but expect the qualitative features to be 
valid for any odd occupancy (see Ref. [H]). We work in 
the experimentally-relevant parameter regime where the 
spin-orbit and orbital-Zeeman couplings are small com- 
pared with the nanotube bandgap and the energy of the 
longitudinal motion in the QD. Here, the longitudinal 
and sublattice orbital degrees of freedom are effectively 
frozen out, leaving behind a nominally four- fold degen- 
erate low-energy subspace associated with the remaining 
spin and valley degrees of freedom (see Refs. [51 I^V 

A simple model describing the spin and valley dynam- 
ics in this low-energy QD subspace, incorporating the 
coupling of electron spin to deflections associated with 
the flexural modes of the CNT [2SJ [5S] , was introduced 
in Ref. |27j . In principle, the deformation-potential 
spin-phonon coupling mechanism [1 11 is also present. 
The deflection coupling mechanism is expected to domi- 
nate at long phonon wavelengths, while the deformation- 
potential coupling should dominate at short wavelengths 
(see discussion in [27^). For simplicity we consider only 
the deflection coupling mechanism, but note that the ap- 
proach can readily be extended to include both effects. 

The Hamiltonian describing this system is [TUl [Ml \T7\ 

H = ^T3{s-t)+AKK'Tl-HorhT3{B-t)+fiBis-B), (1) 



where Ago and Akk' denote the spin-orbit and interval- 
ley couplings, Ti and Si are the Pauli matrices in valley 
and spin space (the pseudospin is frozen out for the states 
localized in a QD), t is the tangent vector along the CNT 
axis, and B denotes the magnetic field. Note that the 
spin-orbit coupling has contributions which are diagonal 
and off-diagonal in sublattice space [TBl - EUl 12^ . When 
projected onto to a single longitudinal mode of the quan- 
tum dot, the effective Hamiltonian given above describes 
the coupling of the spin to the nanotube deflection at the 
location of the dot |24j . 

For a nominally straight CNT we take t pointing along 
the z direction, giving s ■ t = and B • t — Bz- Here 
we find the low-energy spectrum shown in Fig. [2| The 
two boxed regions indicate two different two-level systems 
that can be envisioned as qubit implementations in this 
setup: we define a spin qubit [8] (S) at strong longitudi- 
nal magnetic field, near the value B* of the upper level 
crossing, and a mixed spin-valley or Kramers (K) qubit 
[ID] , which can be operated at low fields applied either in 
the longitudinal (Kz) or perpendicular (Kx) directions. 

We now study how these qubits couple to the quantized 
mechanical motion of the CNT. For simplicity we con- 
sider only a single polarization of flexural motion (along 
the x-direction), assuming that the two-fold degeneracy 
is broken, e.g., by an external electric field. A general- 
ization to two modes is straightforward. 

A generic deformation of the CNT with deflection u{z) 
makes the tangent vector t(z) coordinate-dependent. Ex- 
panding t(z) for small deflections, we rewrite the cou- 
pling terms in Hamiltonian ([T]) as s • t ~ -f [du/dz)sx 
and B • t ~ B^ + {du/dz)Bx. Expressing the deflec- 
tion u{z) in terms of the creation and annihilation op- 
erators and a for a quantized flexural phonon mode, 
u{z) = /(z)^(a -I- a'^), where f{z) and £o are the wave- 
form and zero-point amplitude of the phonon mode, we 
find that each of the three qubit types (S, Kx, Kz) ob- 
tains a coupling to the oscillator mode which we describe 
as 

— = -^a3+ g{a + a^)ai+ujpa^ a + 2\{a + a^) cosut. (2) 

Here the matrices iJi^a are Pauli matrices acting on the 
two-level qubit subspace, and we have included a term de- 
scribing external driving of the oscillator with frequency 
Lo and coupling strength A, which can be achieved by 
coupling to the ac electric field of a nearby antenna [4]. 
Below we describe the dependence of the qubit-oscillator 
coupling g on system parameters for each qubit type (S, 
Kx, or Kz). The derivation of Eq. ^ is detailed in [2i] . 

For the spin qubit (S), the relevant two- fold degree of 
freedom is the spin of the electron itself. Therefore in 
Eq.([2]) we have as = Sz and cti = Sx, and the qubit 
levels are split by the Zeeman energy, measured relative 
to the value B* where the spin-orbit-split levels cross, 
hujq = fisiB — B*). A spin magnetic moment of /is is 
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assumed, and B* « Aso/2/i_B for Akk' ^ Ago- For the 
qubit-resonator coupling, we find g = Aso(/')^o/2'\/2, in- 
dependent of B. Here, (/') is the derivative of the wave- 
form of the phonon mode averaged against the electron 
density profile in the QD. 

For a symmetric QD, positioned at the midpoint of the 
CNT, the coupling matrix element proportional to (/') 
vanishes for the fundamental and all even harmonics (the 
opposite would be true for the deformation-potential cou- 
pling mechanism). The cancellation is avoided for a QD 
positioned away from the symmetry point of the CNT, 
or for coupling to odd harmonics. Here, for concreteness, 
we consider coupling of a symmetric QD to the first vi- 
brational harmonic of the CNT. Using realistic param- 
eter values H m m [50], i = 400 nm, 4 = 2.5 pm, 
Aso = 370 AteV, Akk' = 32.5 ^leY, fi^vh = 1550 ^eV/T, 
and ojp/2'n- = 500 MHz, we find g/27r « 0.56 MHz, irre- 
spective of the magnetic field strength B along the CNT. 

For the Kramers qubits (Kx and Kz), both LUq and 
g depend on B. The qubit splitting for the Kx 
qubit is controlled by the perpendicular field, hujq = 
IJLb{'2Akk' / A)Bx, while for the Kz qubit, it is controlled 
by the longitudinal field hu]q = {^b + A'orb(Aso/A))i32, 
where A = A^^ -I- 4A|^^, denotes the zero-field split- 
ting between the two Kramers pairs. Resonant coupling 
occurs when ujq — iOp. This condition sets the relevant 
value of B^ {B^) in the case of the Kx (Kz) qubit; the 
parameters above yield B^ « 103 mT {Bz « 0.6 mT). 

The qubit-cavity coupling for the Kx qubit increases 
linearly with the applied perpendicular field, Hg — 
-((/')4/V2)(MorbA,o/A + AiBA2jA2)B„ while for the 
Kz qubit it scales with the longitudinal field, hg = 
((/')4/v^)(Morb2AKK'Aso/A2)B^. Using the values of 
Bx and Bz obtained above, we estimate couplings of 
g/2-K « 0.49 MHz for the Kx qubit, and g/2TT w 0.52 kHz 
for the Kz qubit. Thus the coupling for the Kx qubit is 
comparable to that of the spin qubit, while the coupling 
of the Kz qubit is much weaker. Therefore, we restrict 
our considerations to the spin and Kx qubits below. 

Ref . m reports the fabrication of CNT resonators with 
quality factors Q w 150,000. We take Q = 63,000 for 
the following estimate. Together with the oscillator fre- 
quency ujp/2Tr = 500 MHz, this value of Q implies an 
oscillator damping rate of F « 5 • 10^ s^^ ^ g. Because 
of the near-zero density of states of other phonon modes 
at ujq, it is reasonable to assume a very low spontaneous 
qubit relaxation rate 7. These observations suggest that 
the so-called "strong coupling" regime of qubit-oscillator 
interaction, defined as r,7 ^ g, can be reached with 
CNT resonators. 

To quantify the system's response in the anticipated 
parameter regime, we study the coupled qubit-oscillator 
dynamics using a master equation which takes into ac- 
count the finite lifetime of the phonon mode as well as 
the non-zero temperature of the external phonon bath. 



For weak driving, A <C Wp, and ss « w ^ 
we move to a rotating frame and use the rotating wave 
approximation (RWA) to map the Hamiltonian, Eq.([2|, 
into Jaynes-Cummings form [31] 

— — = — cr3+.g(a(T++a'(T_)+Wpa'a + A(a + a'), (3) 

where = uji — lo. Including the nonunitary dynamics 
associated with the phonon-bath coupling, the master 
equation for the qubit-oscillator density matrix p reads: 

i 1 
p = -- [ifRWA,p] + [uB + [apa) - -{a)a,p}) 

+ ubV {a^ pa-]^{aa\p}), (4) 

where ub = Xj (e!^^vl^BT _ j^-j -g -j-j^g bath-mode Bose- 
Einstein occupation factor, and fc^ is the Boltzmann con- 
stant. 

Because of the phonon damping, in the long-time limit 
the system is expected to tend towards a steady state, 
described by the density matrix p. We study these 
steady states, found by setting p = in Eq.Q, using 
both numerical and semiclassical analytical methods. In 
Figs. [3^,c we show the steady-state phonon occupation 
probability distribution P(6uj,n) as a function of the 
drive frequency-phonon frequency detuning duj = —u)p 
and the phonon occupation number n, for the case where 
the qubit and oscillator frequencies are fixed and degen- 
erate, LUq = ujp (see caption for parameter values). Pan- 
els a and c compare the cases with and without qubit- 
oscillator coupling. In Figs.|3jD and[3]i we show the aver- 
aged phonon occupation number n(5u]) — nP{5uj, n), 
which is closely related to the mean squared resonator 
displacement in the steady state: X"^ — x'^ — i'^{n -\- \) . 
For (7 ^ 0, we observe a splitting of the oscillator reso- 
nance, which is characteristic of the coupling to the two- 
level system, and can serve as an experimental signature 
of the qubit-oscillator coupling. For drive frequencies 
near the split peaks, the phonon number distribution is 
bimodal (Fig. |3j') showing peaks at n « and at high-n, 
indicating bistable behavior (see below). 

For strong excitation, where the mean phonon occupa- 
tion is large, we expect a semiclassical approach to cap- 
ture the main features of the system's dynamics [TTJ [35] . 
Extending the approach described in [TT] to include dis- 
tinct values of the qubit, oscillator, and drive frequencies, 
u)q,ujp, and w, we derive semiclassical equations of mo- 
tion for the mean spin and oscillator variables (see [24]). 
The steady-state values of the mean squared oscillator 
amplitude obtained from the resulting nonlinear system 
are shown in Fig. [3J2. In the vicinity of the split peak we 
find two branches of stable steady-state solutions, indica- 
tive of bistable/hysteretic behavior The semiclassical 
results in Fig. |3]3 are in correspondence with the phonon 
number distribution in Fig. |3j:, and explain its bimodal 
character. Similar oscillator instabilities have been used 
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FIG. 3: Response of the spin-oscillator system, (a) Phonon number probability distribution P{n,Suj), (b) average phonon 
occupation n and root mean squared displacement X of the uncoupled driven CNT resonator {g — 0), as functions of the drive 
frequency-oscillator frequency detuning Suj — uj — cUp. The parameters are T — 50 mK, ujp/2-K = 500 MHz, F = 5 • 10^ s"'^ 
and \/2n — 0.027 MHz. The same quantities are plotted in (c) and (d) for a resonantly coupled qubit-oscillator system (i.e., 
ujq = Ljp), with coupling constant g/2n = 0.5 MHz and further parameters as in (a) and (b). (e) Steady-state oscillator 
response from the semiclassical calculation, corresponding to the parameters of (c) and (d). The green solid (purple dashed) 
lines describe stable (unstable) solutions, (f) Bimodal phonon number distribution, taken along the dashed vertical line of (c). 
(g,h) Root mean squared value X of the resonator amplitude in the coupled spin qubit - oscillator system at (g) T = and (h) 
T = 50 mK, as functions of magnetic field detuning SB (detuning the qubit frequency away from resonance with the oscillator) 
and drive frequency-oscillator frequency detuning Slj. 



as the basis for a sensitive readout scheme in supercon- 
ducting qubits [H], and may potentially be useful for 
mass or magnetic field sensing applications where small 
changes of frequency need to be detected. 

To predict the oscillator response to be detected via 
a charge sensor (see below), we solve for the station- 
ary state of Eq. ^ directly for a range of driving fre- 
quencies, qubit-oscillator detunings (set by the magnetic 
field), and temperatures T. In Figs, [sj^g) andjsjjh), we 
show the T — and T = 50 mK root mean squared 
oscillator amplitude X cx y^n + 1/2 as function of mag- 
netic field B and drive frequency, for the case of a spin 
(S) qubit. The value 6B — corresponds to resonant 
coupling Ljq = u)p. These results also apply for the Kx 
qubit, if the magnetic field axis is adjusted appropriately. 
In the zero-temperature case, only half of the eigenstates 
tiw± « hiijp =p hg^ / [ujp — ujq) of Eq. ([s]) can be efficiently 
excited by the drive at fixed SB, giving rise to the upper 
(lower) feature in Fig. [3^' for SB < {6B > 0). However, 
for T > hiOq, both branches of the Jaynes-Cummings lad- 
der can be efficiently excited (Fig.jsji). This is a distinct 
and experimentally accessible signature of the strong cou- 
pling at finite temperature. Note that the vacuum Rabi 
splitting is also observed (see arrows in Fig.jsji), but fea- 
tures arising from nonlinearity in the strongly driven sys- 
tem dominate by more than 2 orders of magnitude. 

Displacement detection of nanomechanical systems is 
possible using charge sensing [51 135] , where the conduc- 



tance of a mesoscopic conductor, such as a QD or quan- 
tum point contact, is modulated via capacitve coupling 
to the charged mechanical resonator. Furthermore, the 
qubit state itself can be read out using spin-detection 
schemes developed for semiconductor QDs or by a 
dispersive readout scheme like that commonly used in 
superconducting qubits coupled to microwave resonators 
[57] . The dispersive regime can be rapidly accessed by, 
e.g., tuning the resonator frequency using dc gate pulses 
which control the tension in the CNT [1]. 

In summary, we predict that strong qubit-resonator 
coupling can be realized in suspended CNT QDs with 
current state-of-the-art devices. The coupling described 
here may find use in sensing applications, and in spin- 
based quantum information processing, where the CNT 
oscillator enables electrical control of the electron spin, 
and, with capacitive couplers, may provide long-range 
interactions between distant electronic qubits [TH] [55] . 
Combined with control of the qubit via electron-spin- 
resonance |39j , the mechanism studied here could be uti- 
lized for ground-state cooling and for generating arbi- 
trary motional quantum states of the oscillator |15| . 
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Note: While completing this manuscript, we became 
aware of a related work [40] that describes the theory of 
the spin-phonon coupling in a CNT resonator QD, and 
its consequences in the spin blockade transport setup. 
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APPENDIX A: EFFECTIVE QUANTUM DOT HAMILTONIAN 

In the following, we derive Hamiltonian (1) of the main text, which describes the spectrum of a quantum dot 
formed in a straight, suspended carbon nanotube, in the absence of phonons (i.e. in a static tube). As usual, we start 
from a tight-binding description of a graphene sheet, which is then rolled up with the condition of periodic boundary 
conditions. Using the conventions as in Weiss et al. [1 , this gives the following Hamiltonian for the longitudinal 
degree of freedom 

Fo = VYP^a2 + A^CTi + 1 • st3(ctiAi + Aq) + V{z). (1) 

Here z and t represent the coordinate and unit vector in the direction of the tube, vp is the Fermi velocity, ai,Ti, and 
Si are Pauli matrices in sublattice, valley and spin spaces, respectively. Note that to translate between the convention 
used here and that used in e.g. Ref. HI cji must be replaced by T3CT1. The energy gap between the valence band and 
the conduction band is 2Ag, where Ag = hvpiy /'^R) + Ac with 2Ac being the curvature induced minigap, which is 
typically of order 10 meV, but proportional to cos 30, where is the chiral angle of the tube. For nominally metallic 
tubes, v ^ 0. The spin-orbit interaction has two terms, one that is diagonal in sublattice (the Aq term) and one 
which is off-diagonal (the Ai term). The spin-orbit interaction connects the spin projection along the tube axis to the 
JC, K' (valley flavor) quantum number through the prefactor (t • s)r3. The microscopic derivation of this Hamiltonian 
can be found in Refs.EHU Finally, the term V{z) describes the confining potential of the quantum dot. It is assumed 
to be smooth on the atomic scale and hence has no sublattice structure. 

Furthermore, an applied magnetic field couples to both the spin and orbital degrees of freedom. The coupling to 
the orbital degree of freedom appears through the Aharonov-Bohm flux, which modifies the boundary condition of 
the circumferential wave vector and hence changes Ag. In total, the Hamiltonian due to magnetic field is 

HB=Horb + Hs, Horb ^ -IJ-B{l/fi)criT3i -B, Hs = ^gfiBS-B, (2) 

where I = mvpR. 

There are several ways to arrive at the Hamiltonian in Eq. (1) in the main text. Here we will assume a hierarchy 
of energy scales, typical of many experimental realizations of nanotube quantum dots, namely 

Ag > Si > , Eso , Ekk' , (3) 

where is the level spacing due to longitudinal quantization and Eb, Ego, and Ekk' are the energy changes due 
to the external magnetic field, spin-orbit coupling, and valley mixing, respectively. This allows us to first solve for 
the dot wavefunction in absence of these three contributions and then project onto a single longitudinal mode. 

In passing we note that in order to get more information about dependence of the orbital magnetic moment and 
spin-orbit coupling on the number of electrons in the quantum dot, one has to be more precise and use a specific form 
of the confining potential, e.g. assuming a square well potential, as in Refs. [TJ [71 This was done in Refs. [SI [51 where 
a method to experimentally extract the two spin-orbit parameters Aq and Ai, as well as /Zorb, was shown. 

Now imagine that one has solved for the case without magnetic field, spin-orbit coupling, and mixing between K 
and K' . This gives a set of longitudinal wavefunctions, each one four-fold degenerate due to the spin and valley 
degrees of freedom. The energy splitting between these shells is El ^ -Es, i?sO: Ekk' ■ We label these states by the 
valley and spin quantum numbers r = ±1 and s = ±1, respectively, which indicate corresponding eigenvalues under 
T3 and Sz- Projected onto eigenstates of the spatial coordinates z and c (the circumferential coordinate), the wave 
function in the envelope-function representation has the form[T] 



{{z, c)|t, s) = — =.0(z)?7<, ®Xs® Xr, 



(4) 
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where k± = is the wave vector associated with the gap, r]cr = l)""" is a pseudospinor describing the 

sublattice degrees of freedom, Xs and Xr describe the spin and valley degrees of freedom, respectively, x+ = (1 O)"*" 
and X- — (0 ^)'^ i and the precise form of the envelope wave function (j>{z) depends on the confining potential. 

We can now take matrix elements with respect to Hb and i?so = t • 8X3(0' j^Aj^ + Aq) and also a term describing the 
microscopic disorder that couples valleys: Hkk' = Vkk' , where Vkk' is a short-range disorder potential depending 
on the longitudinal and circumferential coordinate operators. (Note that short-range disorder can be systematically 
incorporated into the envelope-function description, see e.g., Refs. 13 [9l) This procedure will produce a Hamiltonian 
of the form in Eq. (1) in the main text, 

H = ^■^3(8 • t) + Aa'A"T"1 - AiorbT3(B • t) + iib{s ■ B), 

with 

A,o = 2(+|f7iAi -f A0I+), Akk' = {-\Vkk'\+), Morb = ^^(+ki|+). (5) 

Here the kets |±) stand for the orbital states with r = ±1. Note that in general, the valley-mixing term can include 
both Ti and T2, but an appropriate unitary transformation in the valley space can be used to put it into the form 
above with real Akk'- 



A.l Derivation of the coupling to vibrations 

Next we look at how the vibrations couple to the four states of the quantum dot. The amplitude of the tube is in 
terms of the harmonic oscillator raising/lowering operators given by 

«(z) = /(z)^(a + at), (6) 

where we focus on a single vibrational mode. The coupling to the spin is via the change of the tangent direction given 
by 

6t -—X, (7) 
dz 

where x is perpendicular to the tube and in the plane of the vibration. The interaction Hamiltonian then becomes 

-ffs,vib = (5t • S T^{aiA^ + Aq) - /iorbTgCTlfJt • B. (8) 

As above, we now project onto a single longitudinal mode, thus taking matrix elements of Hsyih in the basis |t, s). 
Such matrix elements involve form factors like 

{T,s\aJ'{z)\T\s') =Srr'{r,s\aJ'iz)\T,s') =Srr'Sss'F,,r. (9) 

At this point we note that the coupling is small for a symmetric dots and even harmonics, because the F factors then 
tend to cancel, see discussion in the main text. The effective Hamiltonian for coupling between the four states of the 
quantum dot and the vibration now becomes 

-ffs.vib = ^ (a + a^) {s^T^iFiA^ + FqAq) ~ B^r^^iorhFi} . (10) 

We see from Eq. ( [lO| ) that the coupling of the vibrations to quantum dot states have different form factors from 
what one would get by simply setting Q into Eq. (1) of the main paper. However, when the energy scales are clearly 
separated as in (|3| the eigenstates |r, s) are eigenstates of ai (which can be seen from ([!])) and therefore Fq — Fi 



(for the conduction band). Therefore, we do not need to take into account the different form factors in (10), which 
simplifies the analysis and we can write Fq = Fi = (/'). In this language Eq. (|lO|) becomes 



HsMh = ^ + {^^jT-sAso - SxTsAiorb} (/') , (H) 



which is the result used in the main text. 



As mentioned the expression (11) was derived under the assumption that the gap dominates over the longitudinal 
size quantization energy, which is valid for few-electron quantum dot. However, it is important to note that one could 
easily extend this to the more general case at higher energies by including the difference in form factors Fq and Fi 
without changing the conclusions and structure of our results qualitatively. 
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APPENDIX B: QUBIT-PHONON COUPLINGS 

We treat three different qubit realizations in the main text: the spin qubit (S), the Kramers qubit in a magnetic 
field perpendicular to the carbon nanotube (CNT) (Kx), and the Kramers qubit in a parallel-to-CNT magnetic field 
(Kz). Below, we express the three qubit Hamiltonians Hg, Hkx and Hkz as functions of system parameters and CNT 
deformation. To clarify the correspondence with the Jaynes-Cummings Hamiltonian in Eq. (2) of the main text, we 
list the formulas for the qubit frequency, the qubit-phonon coupling, as well as numerical estimates for the latter, in 
Table H 



B.l Spin-phonon coupling 

At the finite value 




of a longitudinally-applied magnetic field, the quantum dot (QD) energy spectrum shows a crossing of the energies of 
a pair of spin states belonging to the same valley (see Fig. 2 of the main text). Around this point these two levels 
are energetically well separated from the other QD levels. We call this two- level system the spin qubit (S). If the 
dynamics is restricted to these two levels, it can be described by the following effective Hamiltonian: 

=Mi3S.(S-S*) + ^^s, (13) 

where we assume that the effect of valley-mixing is negligible, ^kk' = 0- Averaging over the z coordinate using the 
charge density n[z) of the electron occupying the CNT QD yields 

Hs^ J dzniz)H,{z)^^iBsAB-B*) + ^s,J dzn{z)'^^^ ^iBS.{B ~ B*) + 1^^^ . (14) 



B.2 Kramers qubit-phonon coupling in a perpendicular magnetic field 

At zero magnetic field, the ground state of the CNT QD, i.e., the ground state of the Hamiltonian H in Eq. (1) 
of the main text, is formed by a pair of time-reversed states (Kramers pair). The twofold degeneracy is maintained 
even in the presence of spin-orbit interaction and valley mixing. At small enough magnetic field these two states split 
up, but they remain energetically well separated from higher-lying states. We call this two-level system the 'Kramers 
qubit' [in] in a perpendicular field (Kx). [Similar considerations hold for the first excited Kramers pair, i.e., the two 
higher-lying energy eigenstates of H in Eq. (1) of the main text.] These energetically split states, in the absence of 
CNT deformation, will be denoted here as |-|-) and |— ). Starting from the Hamiltonian H in Eq. (1) of the main text, 
averaging over z using the electron density n(z), and incorporating the effect of the two higher- lying states on |-|-) and 
|— ) via a second-order Schrieffer- Wolff transformation, we find that the dynamics restricted to the Kramers qubit in 
the presence of an external magnetic field B — B^^i. and CNT deformation u{z) is described by the Hamiltonian 



Hkx — Bx 



CT3- 



2mbA 



KK' 




S 



^'-^KK' 



Here, (Ti_3 are the Pauli matrices in the qubit basis, i.e., 0-3 = |+)(+| — |— )(— | and ai = |+)(— | + |— )(+|- 



(15) 



B.3 Kramers qubit-phonon coupling in longitudinal magnetic field 

In the absence of CNT deformation, a parallel-to-CNT magnetic field splits two low-energy Kramers doublet of 
H in Eq. (1) of the main text. We call this two-level system the Kz qubit, and denote the two qubit states as |-|-) 
and |— ) in this subsection. Starting from the complete Hamiltonian H in Eq. (1) of the main text, averaging over z 
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Hamiltonian 




kg 


(7/27r (numerical) 


Hkx 
Hkz 




I1b{B-B*) 






( MorbAso 1 A'fjA^^ ^ 




U.tJU iVlllZi 

0.49 MHz 
0.52 kHz 




B 


W^-+4A^^, A^„+4A^^^J 

^{f')ta Morb2A^^,A,o 
Ai„+4A^^^, 



TABLE I: Correspondence between the terms of the Hamiltonian in Eq. (2) of the main text and those of the qubit-phonon 
Hamiltonians derived for the three different qubits. The parameters used to calculate the last column are L — 400 nm, £0 = 2.5 
pm, Ago ~ 370 /ieV, Akk' ~ 32.5 pieV, /lorb ~ 1550 fieV/T. The estimate {/') = 2^/2/ L has been used (see text). 



using the electron density n{z), and applying a second-order SchriefFer- Wolff transformation to describe the effect of 
the higher-lying Kramers pair to the Kz qubit, we find that the dynamics of the latter in the presence of an external 
magnetic field B = B^z and CNT deformation u{z) is described by the Hamiltonian 



H 



Kz 



/ du\ ^orbSA^^'A 

(7l( - 



\dz/ AL + 4A2 



/iorb^s 



■ 4A2 



(16) 



As before, CTi^s are the Pauli matrices in the qubit basis, i.e., (T3 — |+)(+| — |— )(— | and CTi — |+)(— | + |— )(+|- 



B.4 Deformation of the CNT 



All three qubit-phonon Hamiltonians Hg, Hkx and Hkz resemble the Jaynes-Cummings Hamiltonian of cavity 
quantum electrodynamics. This becomes more apparent if we express the z-dependent displacement u{z) in terms of 
phonon annihilation a and creation operators: 



uiz) = /(z)-^(a-f a^). 



(17) 



Here f(z) is a dimensionless function describing the shape of the standing- wave bending phonon mode under consid- 
eration (normalization: J p{z)dz = L) and £0 is the ground-state displacement of that mode. 

As apparent from Table ll the qubit-phonon coupling vanishes if (? oc (/') = J dz^^ljp-n{z) = 0. Therefore, g vanishes 
if the setup is perfectly left-right symmetric along the CNT (z) axis and a bending mode with an even number of 
nodes is considered. To have a finite qubit-phonon coupling, either the left-right symmetry of the setup must be 
broken or a flexural mode with odd number of nodes should be considered. In the main text and also here we treat 
the second case: we investigate the coupling of the first excited flexural phonon (1 node) to the various qubits. 

The displacement field of the first harmonic of the resonator can be approximated by 



fiz) 



V2sir 



2tt 

T 



z + 



L 



(18) 



where we assume that the CNT is suspended at points z = —L/2 and z = L/2. Approximating the charge density 
with a step function symmetrically covering the length fraction ^ of the suspended part of the CNT, we obtain 



2tt 



z + 



2^2 sin TT^ 



The fraction sin(7rf)/^ is ~ 1 if ^ is smaller than 1, and therefore we make the approximation (/') 
numerical estimates appearing in the main text and in Table |Tj 



(19) 

2V2/L in the 



APPENDIX C: SEMICLASSICAL EQUATIONS OF MOTION 



In this section we develop semiclassical equations of motion for the coupled qubit-oscillator system, valid in the 
regime of large oscillator excitation. We follow the procedure of ReffTTl this time allowing for different qubit (wg). 
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oscillator (wp), and drive (oj) frequencies. Our aim will be to find a closed set of equations for the time dependence of 
the expectation values of the oscillator and qubit coordinates, (a), (cr~), and (cr^). We evaluate the time derivatives 
of these observables using ^ (O) = Tr[p0], with the time-dependence of the density matrix given by Eq.(4) of the 
main text (below we set h = 1), 

P = ~i[HBWA,p\ + (ns + l)r (apa^ - \{a^a,p}] + ubT [ pa - \{aa\ p} 



with ub = l/(e"'''^ — 1) and i?RWA = + 9{a(J+ + a' cr_) + uipa'a + A(a + a^). Here, in the rotating frame, the 
reduced frequencies are given by uJi = uJi — lj. The qubit-phonon coupling is denoted by g, and the strength of the 
external driving field is denoted by A. 

Using the commutation rule aa^ = a^a + 1 and the cyclic property of the trace, Tr[74_B] — Tr[_Byl], we find: 

(a) = {—iCjp — T/2){a)—i\ — ig(^a^^ 
(<t") = -iujq (a^) + ig (aa^) 
i&s) = -2tg{{aa+)-{a^a-)). 

We close this set of equations by neglecting correlated fluctuations between the qubit and oscillator degrees of freedom. 



factoring the averages as (acrs) 



and 



{cr~^). Using (cr"'") = (cr )* and (a'^) = (a)*, and as in 



Ref llll defining the complex variables z = (a) and v = 2 {a ), and a real variable m = (trs), we obtain 



V 

m 



- {iujp 4 

-iujqV - 
-ig{zv^ 



r/2)z- 

- 2igmz 



iX 



vz* 



(20) 

(21) 
(22) 



In this representation, the real and imaginary parts of z describe the oscillator coordinate and momentum, the complex 
variable v describes the x and y Bloch vector components of the qubit state, and m describes the qubit polarization. 



The dynamics described by the nonlinear system in Eqs. (|20|)-(|22|) can be quite complex. Here we focus on steady 
state solutions, z = 0,v 
values, we obtain 



0,m = 0. Setting w = in Eq. (21), and introducing over-bars to indicate steady state 

- 2g__ 



(23) 



Note that m = is automatically satisfied under this condition, see Eq. (22 1 



Within the semiclassical description, and in the absence of decoherence acting directly on the spin, the variables v 

, , _ / . 2 _N -1 

and m describe a vector of unit length, -|- m 
the square root of both sides gives 



1. Using Eq. (23 1 for we find m 



1 



Taking 



m± = ± 1 



-1/2 



(24) 



where the subscript ± indicates two branches of solutions to the square root. Setting i = in Eq. (20), and using 



Eqs. (23) and (24), we find that the oscillator amplitude in the steady state satisfies the relation 

2 



A2 



I2±l 



= (r/2)^ 



UJp ± 



(25) 



Self-consistent solutions to Eq. (25 1 can easily be found numerically. In many parameter regimes, multiple solutions 



exist due to the non-linearity introduced by the qubit-oscillator coupling. As shown in Fig.3e of the main text, the 
branches of stable fixed points match well with the peaks in the phonon number distribution, indicating the utility of 
the semi-classical approach. 

Along with the presence of multiple solutions, we expected the typical manifestations of multistable behavior, such 
as hysteresis and sharp instabilities. The sensitivity to the steady state oscillator amplitude near instability points 
where stable steady-state solutions disappear may be useful for sensing applications. Closely related behavior has 
already proved quite useful in providing a sensitive read-out mechanism for superconducting aubits.[12j 
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APPENDIX D: INTERPRETATION OF RESULTS FOR CHARGE-SENSING-BASED DETECTION 



As stated in the main text, the oscihatory motion of the CNT resonator can be detected using a charge sensing 
scheme in which the conductance of a mesoscopic conductor, such as a QD or quantum point contact, is modulated via 
capacitive couphng to the charged mechanical resonator. At a given source-drain bias on the mesoscopic conductor, 
the current depends on the displacement u of the resonator due to their capacitive coupling. A nonlinear dependence 
of the current on the resonator displacement is desired for displacement sensing, I{t) sa /o + iMt) + hu^it). Such 
a dependence is present in, e.g., a QD tuned to the middle of a Coulomb-blockade peak (corresponds to /i = and 
I2 < 0), or to the onset of a Coulomb-blockade peak [Ii > and /2 > 0). Here we show that under such conditions, 
and neglecting detector back-action on the oscillator, the steady-state time-averaged current (/) through the charge- 
sensing mesoscopic conductor is sensitive to the steady-state average number n of phonons in the oscillator. Therefore, 
the results plotted in Fig. 3b, d, g, h of the main text can be interpreted as being proportional to the measured signal 
(/) — /o in the charge-sensing setup described above, hence that setup would allow for the experimental confirmation 
of the predicted features. 

The steady-state time-averaged current (/) through the charge-sensing conductor is given by 



(/) = lim - (dtl[t) = lo+h lim ^ /di 
T^oo r Jo T Jo 



dtv?{t). 



(26) 



The instantaneous square displacement is expressed with the steady-state density matrix p and the position operator 
(i), both represented in the rotating frame, as 



X 



u^{t) = Tr{pa;2(i)} = Tr <^ p 



V2 



(27) 



Substitution of this expression to Eq. ( 26 1 yields 



(/) - Jo = h^l 



(28) 



as stated above. 
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